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Abstract: VERITAS is an array of four 12-m diameter imaging atmospheric-Cherenkov telescopes located in southern 
Arizona. Its aim is to study the very high energy (VHE: E > 100 GeV) 7-ray emission from astrophysical objects. In 
addition to the study of blazars, the VERITAS extragalactic science program develops a comprehensive observational 
program of extragalactic non-blazar sources. The study of Active Galactic Nuclei (AGN) is intensively pursued through 
the large MWL observational campaign on non-blazar radio galaxies. The success of these MWL campaigns has for the 
first time provided insights to the inner structures of jets responsible for 7-ray emission. The problem of the origin and 
acceleration of ultra-high energy cosmic rays is pursued through the observation of the starburst galaxy M 82, whose 
detection at VHE delivers important insights on the possible acceleration mechanisms. Finally, investigation of globular 
clusters places important limits on the millisecond pulsar contribution to their addressed 7-ray emission. The VERITAS 
extragalactic non-blazar science program and its results are presented. 
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1 Introduction 

The largest population of VHE-detected 7-ray sources 
are blazars, a sub-category of AGNs in which the ultra- 
relativistic jet, produced by the accretion of matter around 
a super-massive black hole, is aligned within a few de- 
grees to the observer's line of sight Q). The most com- 
monly accepted model to account for the 7-ray emission 
from the jet of AGNs is inverse-Compton scattering of the 
synchrotron photons produced by shock-accelerated elec- 
trons and positrons within the jet itself [2]. In these aligned 
sources relativistic beaming substantially boosts the ap- 
parent flux. Non-blazar AGNs are typically oriented at 
larger angles from the observer's line of sight, becoming 
much more challenging. However, the misalignment en- 
ables imaging of the jet's structure, crucial to identifying 
the emission regions and probing models of the accelera- 
tion mechanism. Given the typical angular resolution of 
the order of several arcminutes in 7-ray instruments, jet 
substructures are not resolved in the 7-ray energy band, but 
they are in other wavelengths. Correlation studies through 
coordinated multi-wavelength (MWL) observational cam- 
paigns on radio galaxies are a viable strategy to investigate 
the physical processes at work in the substructures of the 
jet. 

The advantage of observing radio galaxies is that it is pos- 
sible to study also the rich environment in which they are 
typically located. It has been seen that radio galaxies are 
preferentially located in cluster of galaxies Q. Their pow- 



erful jets energize the intra-cluster medium through the ter- 
mination shocks accompanied by particle acceleration and 
magnetic field amplification. Large scale AGN jets and 
cluster of galaxies are believed to be potential accelerator 
for cosmic rays |4|, therefore the modeling of the dynamics 
of both populations is of particular interest for the cosmic- 
ray community. 

Beside AGN-related environments, starburst galaxies are 
also good candidates as ultra-high energy cosmic rays ac- 
celerators. The active regions of starburst galaxies have a 
star formation rate about 10 times larger than the rate in 
normal galaxies of similar mass, with a consequent higher 
rate of novae and supernovae. The cosmic rays produced in 
the formation, life, and death of their massive stars are ex- 
pected to eventually produce diffuse gamma-ray emission 
via their interactions with interstellar gas and radiation. 

Finally, globular clusters are the closest extragalactic struc- 
tures whose physics is interesting to the 7-ray community. 
They can host hundreds of millisecond pulsars which can 
accelerate leptons at the shock waves originating in colli- 
sions of the pulsar winds and/or inside the pulsar magne- 
tospheres. Energetic leptons diffuse gradually through the 
globular cluster. Comptonization of stellar and microwave 
background radiation is therefore expected to be responsi- 
ble of 7-ray emission. 

The indirect search for dark matter (DM) candidates, is also 
part of the VERITAS extragalactic non-blazar program. 
A dedicated contribution on the VERITAS DM program 
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is presented in a separate proceeding. Highlights on the 
research topics and results of the VERITAS extragalactic 
non-blazar science program are here presented. 

2 The VERITAS Instrument 

The VERITAS detector is an array of four 12-m diame- 
ter imaging atmospheric Cherenkov telescopes located in 
southern Arizona 0. Designed to detect emission from 
astrophysical objects in the energy range from 100 GeV to 
greater than 30 TeV, VERITAS has an energy resolution 
of ~15% and an angular resolution (68% containment) of 
~0.1° per event at 1 TeV. A source with a flux of 1% of the 
Crab Nebula flux is detected in —25 hours of observations, 
while a 5% Crab Nebula flux source is detected in less than 

2 hours. The field of view of the VERITAS telescopes is 
3.5°. For more details on the VERITAS instrument and the 
imaging atmospheric-Cherenkov technique, see ■ 

3 The Extragalactic non-Blazar Science Pro- 
gram 

3.1 Radio Galaxies 

Most of the VERITAS observations of radio galaxies are on 
the radio galaxy M 87. This AGN is located in the center 
of the Virgo cluster at a distance of — 16 Mpc and is cur- 
rently the brightest detected VHE radio galaxy. M 87 was 
originally detected with marginal significance by HEGRA 
at TeV energies Q, and later also by HESS @, VERI- 
TAS |9) and MAGIC [13). This giant radio galaxy has al- 
ways been of particular interest because its jet lies at -20° 
respect to the line of sight and its core and the structure of 
the jet are spatially resolved in X-ray, optical and radio ob- 
servations, thus it is an ideal candidate for correlated MWL 
studies lfl6l . 

In 2008 VERITAS coordinated an observational campaign 
with two other major VHE observatories (MAGIC, HESS), 
overlapping with VLBA radio observations 02). Three 
Chandra X-ray pointed observations have also been per- 
formed during the first half of 2008. Multiple flares at VHE 
have been detected. In X-rays, the inner-most knot in the 
jet (HST-1) was found in low state, while the core region 
was in high state since 2000. Progressive brightening of 
the core region in radio was also seen along the VHE flare 
development. This is an indication that the 7-ray emission 
originates from a region close to the core rather than from 
more distant regions. 

In April 2010, during the seasonal monitoring of M 87, 
VERITAS detected another flare with peak flux of —20% 
of the Crab Nebula flux. During the six-month observation 
period, M 87 was detected at a level of 25.6a above the 
background, with an average flux above 350 GeV equiv- 
alent to 5% of the Crab Nebula flux. Dedicated analy- 
sis in 20-minute bins has been performed on the April 
2010 flaring episode. A spectral analysis has been done 
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Figure 1: (upper plot) VERITAS light curve of the 2010 
seasonal monitoring campaign, (lower plot) Spectral anal- 
ysis for three phases of the April 2010 flaring event: rising 
phase ( circles), peak (squares) and decreasing phase (tri- 
angles). 



on three different phases of the flaring episode: the ris- 
ing phase, the peak and the falling phase. A power-law 
fit has been applied to each phase, showing a hint of spec- 
tral variability: T lise = 2.60 ± 0.31, r pcak = 2.19 ± 0.07, 
r faU = 2.62 ±0.18. 

Other radio galaxies observed by VERITAS include 3C 111 
and NGC 1275. A preliminary analysis of 1 1 hr of quality- 
selected data of 3C 111 results in a flux upper limit of 
- 3% Crab flux above 200 GeV. NGC 1275 is an un- 
usual early-type galaxy located in the center of the Perseus 
cluster. Its radio emission is core dominated, but emission 
lines are also seen, making it difficult to classify it accord- 
ing to the Faranhoff & Riley (FR) classification iflOll . In 
Fall 2008 the Fermi 7-ray space telescope reported the de- 
tection of 7-ray emission from a position consistent with 
the core of NGC 1275. VERITAS observed the core re- 
gion of NGC 1275 for about 11 hr between 2009 January 
15 and February 26, resulting in 7.8 hr of quality-selected 
live time. No 7-ray emission is detected above the analy- 
sis energy threshold of — 190 GeV, resulting in a flux upper 
limit incompatible with the extrapolation of the Fermi-LAT 
spectrum. Under the assumption of a SSC emission mech- 
anism, the VERITAS result suggests the presence of a cut- 
off in the sub- VHE energy range ATI . The detection in 
Summer 2010 of VHE 7-ray emission by MAGIC Q2I has 
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eventually included NGC 1275 among the few interesting 
radio galaxies for future VHE investigation. 

3.2 Clusters of Galaxies 

Observation of clusters of galaxies is done unavoidly dur- 
ing the observation of many radio galaxies. This is the case 
for NGC 1275 and M 87 where the Perseus and Virgo clus- 
ters respectively are observed during the radio galaxy ob- 
servation. However, up to now a dedicated study of the 
clusters themselves has been done only on the Coma clus- 
ter. The Coma cluster is a nearby cluster of galaxies which 
is well studied at all wavelengths |14|. It is at a distance 
of 100 Mpc (z = 0.023) and has a mass of 2 x 10 15 M©. 
Its X-ray and radio features suggest the presence of accel- 
erated electrons in the intergalactic medium emitting non- 
thermal radiation. Beside relativistic electrons, there may 
also be a population of highly energetic protons. Both 
high energy electrons and protons are known to be able 
to produce VHE photons. A total of 19 hr of data have 
been recorded between March and May 2008. No evidence 
for point-source emission was observed within the field of 
view and a preliminary upper limit of ^3% of the Crab flux 
is given for a moderately extended region centered on the 
core 1 15 1. 

3.3 Starburst Galaxies 

M 82 a prototype small starburst galaxy, located approxi- 
mately 3.7 Mpc from Earth, in the direction of the Ursa Ma- 
jor constellation. M 82 is gravitationally interacting with 
its nearby companion M 8 1 . This interaction has deformed 
M 82 in such a way that an active starburst region in its 
center with a diameter of ^1000 light years has been de- 
veloped [19 20 1 . Throughout this compact region stars are 
being formed at a rate approximately 10 times faster than 
in entire "normal" galaxies like the Milky Way. Hence the 
supernovae rate is 0. 1 to 0.3 per year [21 22 1. The high star 
formation rate in M 82 implies the presence of numerous 
shock waves in supernova remnants and around massive 
young stars. Similar shock waves are known to accelerate 
electrons to very high energies, and possibly ions too. The 
intense radio-synchrotron emission observed in the central 
region of M 82 suggests a very high cosmic-ray energy 
density, about two orders of magnitude higher than in the 
Milky Way (23). Acceleration and propagation of cosmic 
rays in the starburst core are thus expected to be responsible 
for VHE 7-ray emission. Theoretical predictions include 
significant contributions from both leptonic and hadronic 
particle interactions. Cosmic -ray ions create VHE gamma 
rays through collisions with interstellar matter, producing 
7T° which decay into 7-ray s. Alternatively, accelerated 
cosmic-ray electrons may inverse-Compton scatter ambi- 
ent X-ray photons up to the VHE range 11201 l24l l25l l26l . 

VERITAS observed M 82 for a total of 137 hours of 
quality-selected live time between January 2008 and April 
2009 at a mean zenith angle of 39°. An excess of 91 
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Figure 2: VERITAS 7-ray flux of M 82 compared to the- 
oretical predictions. The data are given by open diamonds 
with la statistical error bars, and is fitted (x 2 /ndf = 0.1) 
with a power-law function (thick gray line), dJV/d-E oc 
(£/lTeV)- r , with r = 2.5 ± 0.6 stat ± 0.2 syst . The VER- 
ITAS flux upper limit (99% confidence level ll29l ) shown 
at 6.6 TeV is above the extrapolation of the fitted power- 
law function at these energies. The thin lines represent a 
recent model ll25l for the 7-ray emission from M 82. The 
thin solid line is the total emission predicted. The dashed 
lines represent components from ir° decay, and from radi- 
ation from cosmic -ray electrons through IC scattering and 
Bremsstrahlung. The markedly different spectral slopes of 
these dominant components should be noted. 

gamma-ray-like events (~0.7 photons per hour) are de- 
tected for a total 4.8<r statistical post-trials significance 
above 700 GeV. The observed differential VHE gamma- 
ray spectrum is best fitted using a power-law function with 
a photon index V = 2.5 ± 0.6 s t a t ± 0.2 sys (fig©. Compar- 
ison of the VERITAS VHE spectrum with predictions of 
the theoretical models supports a hadronic scenario as the 
dominant process responsible for the VHE emission (27). 

3.4 Globular Clusters 

M 15 (aka NGC 7078) is a very compact (R c = 0.19 pc) 
globular cluster located at a distance of 9.4 kpc. It be- 
longs to the class of "core-collapsed" globular clusters and 
is thought to contain a ~ 2000 M black hole in its cen- 
ter. M 15 contains 7 known millisecond pulsars and a low- 
mass X-ray binary (LMXB) system. That are expected to 
emit 7-rays. Pulsars produce relativistic strongly magne- 
tized winds which, when colliding with other winds, create 
relativistic shocks. Models presented in [28 1 consider the 
case of relativistic leptons injected by the pulsar in such 
relativistic shocks with a power-law energy distribution. 

VERITAS observed the globular cluster M 15 between 
June 10 and June 21, 2010 for a total of 6.2 hours of good- 
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Figure 3: VERITAS 95% c.l. upper limit on M 15. The 
three power-laws represent the the three different funtions 
used to calculate the decorrelation energy. It can be seen 
that the dependency of the flux upper limit on the as- 
sumed spectral index is minimal at the decorrelation en- 
ergy. Models representing different leptonic spectra [28| 
are also shown: 30 TeV cutoff (thick lines); 3 TeV cutoff 
(thin lines); T — —2.1, E m i n = 100 GeV (solid lines); 
r = —2.1, E min = 1 GeV (long-dashed lines); T = —3, 
E m i n = 100 GeV (long-dashed-dotted lines); T — —3, 
E m i n = 1 GeV (short-dashed lines). 

quality live time. The data have been processed with the 
standard analysis procedure and cuts optimized for a weak 
Crab-like source (index T = —2.5, 3% Crab Nebula flux). 
No signal is detected at the cluster's core coordinates. A 
95% confidence level upper limit using the Rolke algo- 
rithm 1 29 1 is calculated assuming a Gaussian background. 
The differential limit is calculated at the decorrelation en- 
ergy, i.e. the energy that minimizes the dependency of the 
flux upper limit on the assumed spectral index. Figure [3] 
shows the VERITAS upper limit superimposed on the mod- 
els explained in [28 1. Three functions used to calculate the 
decorrelation energy are shown as well. It can be seen that 
the 95% confidence level VERITAS upper limit places con- 
straints on the possible scenarios for the spectral energy 
distribution of the accelerated leptons in the pulsar winds. 

4 Conclusions 

The VERITAS extragalactic non-blazar science program is 
well established. In addition to complementing the blazar 
program on the study of the AGN physics, non-AGN 
sources are also investigated. A dedicated VHE study 
on the Coma cluster of galaxies resulted in a flux upper 
limit on the extended region centered on the core. The 
first detection of 7-ray emission from a starburst galaxy 
established a connection between cosmic -ray acceleration 
and star formation. Observation of globular clusters 



provides constraints on predictions of 7-ray emission from 
ms pulsars. 
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